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Specification  of  the  Thickness  of  the 
Topside  of  the  Ionosphere 


1.  INTRODUCTION 

The  electron  density  of  the  upper  ionosphere  is  usually  modeled  by  empirical 
profiles  that  decrease  monotonically  with  height.  The  calculated  values  of  elec- 
tron density  are  critically  dependent  on  the  initial  choice  of  a thickness  parameter 
for  the  interval  just  above  the  peak  of  the  F2  region.  Since  there  are  so  few  field 
observations,  it  is  common  practice  to  equate  the  thickness  of  the  top  and  the  bot- 
tom of  the  F2  region  and  then  to  refer  to  climatological  sources,  such  as  those 

1 2 

given  by  Barghausen  et  al  or  CCIR  to  obtain  an  estimate  of  a topside  thickness 
parameter.  Unfortunately,  few  climatological  descriptions  of  the  thickness  of  the 
lower  ionosphere  are  based  on  detailed  electron  density  profiles. 

In  this  paper,  another  approach  is  suggested:  To  obtain  an  estimate  of  the 
height  of  the  maximum  of  the  F2  region  and  then  to  use  it  to  specify  a thickness 
parameter  for  the  topside.  This  implies  using  some  model  relationship  such  as 
the  one  presented  here. 


1.  Barghausen,  A.  F. , Finney,  J.W. , Proctor,  L.  L. , and  Schultz,  L.  D.  (1969) 

Prediction  Long-Term  Operational  Parameters  of  High-Frequency  Sky-Wave 
Teiecommimlcwions  Systems,  iSisA  Technical  Report  Lku  llO-lTS-7a. 

2.  The  International  Radio  Consultative  Committee  (CCIR)  (1970)  Report  No. 

252-2,  Documents  of  the  Xllth  Plenary  Assembly,  New  Delhi  (ITV,  Geneva). 


As  a pilot  study,  a simple  specification  model  has  been  constructed  from 
profiles  of  the  F region  (1968-1971)  kindly  furnished  by  John  Fvans  of  Lincoln 
[laboratory.  Results  indicate  that  the  mean  topside  thickness  parameter  can  be 
estimated  to  within  about  15  percent,  assuming  that  the  height  of  the  F region  is 
known,  and  that  a useful  reduction  can  be  made  in  the  day  to  day  variability. 


2.  DF.VELOPMKNT  OF  A SIMPLIFIED  MODEL 

3 

First  define  Q as  the  quarter  thickness  of  a parabola  fit  to  the  underside  of 

the  F2  region  and  define  h^^^^  as  the  height  of  the  F2  maximum.  Electron  density 

profiles  resulting  from  the  reduction  of  bottomside  vertical  incidence  ionograms 

4 5 6 

such  as  those  reported  by  Beckei}  Clarke  and  Hammond,  and  Wright  can  be  used 

to  show  that  the  mean  noon  and  mean  midnight  values  of  vary  with  season  and 

solar  cycle  in  the  same  sense  as  the  mean  value  of  h (Figure  1). 

max 

Experimental  profiles  of  the  region  above  h are  few  and  scattered.  Avail - 

max 

able  for  this  study  were  the  archive  profiles  obtained  from  Evans  which  were  pro- 
duced by  reduction  of  the  incoherent  scatter  profiles  from  the  Millstone  Hill  radar. 
About  151  days  spanned  the  period  July  1968  through  December  1971  (about  3400 
profiles  edited). 

A parabola  was  fitted  both  to  the  top  and  bottom  of  the  individual  profiles.  Let 
usdefineQ^  as  the  parabolic  quarter  thickness  for  the  topside  segment  between 
and  a higher  height  where  the  electron  density  equaled  0. 7 of  the  density  at  the  F2 
peak.  Similarly  define  Q for  the  bottomside  segment  between  h and  either  the 
height  whose  density  is  0.7  of  the  peak  density  or  200  km,  whichever  is  greater. 

The  bchavic  ! f vs  is  illustrated  in  Figure  2 for  thi’ee  time  samples; 

day,  night,  and  a transition  period  delimited  by  plus  and  minus  two  hours  about 
ground  sunrise  and  ground  sunset.  All  observations  from  the  15  1 days  are  included. 
Count  is  given  in  each  cell  where  size  is  2 km  in  quarter  thickness  by  5 km  in  h 

max 

The  agreement  of  the  variation  of  quarter  thickness  with  h in  Figures  1 and  2 

m cix 

that  an  empirical  model  derived  from  the  limited  sample  of  incolierent 
scatter  data  might  be  useful  over  the  wider  latitude  and  time  range  of  the  data  of 
Figure  1.  If  the  data  are  considered  in  small  blocks  delimited  by  botti  time  and 
season,  tlien  a linear  relation  such  as 

PigROtt,  W.  R. , and  Rawer,  K,  (1972)  URSI  Handbook  of  lonogram  Interpreta- 
tion  and  Reduction,  2nd  ed. , Report 'VAG-23.  ^ ^ 

4.  Becker,  W.  (1970)  The  standard  profile  of  the  mid-latitude  F region  of  the 

^nosphere  as  deduced  from  bottomside  and  topside  ionograms.  Space 
Research  XII.  Akademie  Verlag,  Berlin,  1241-1252 

5.  Clarke,  C. , and  Hammond.  E.  (1965)  J.  Atmos  Terr.  Phys.  27:551. 

6.  Wright,  J.W.  (1962)  J.  Research  Natl.  Bur.  Standards  66D;297. 
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Figure  1.  Parabolic  CJuarter  Thickness  of  Bottom  of  F2 
Region  vs  Height  of  Peak  of  F2  Region.  Seasonal  medians 
for  several  latitudes  and  for  solar  cycle  extremes 
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Fiuure  2a.  Parabolic- 
Quarter  Thickness  of  Top 
of  F2  Region  vs  Height  of 
Peak  of  F2  Region.  Individ- 
ual observations  from  incoherent 
scatter  at  Millstone  Hill:  daytime 
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Figure  2b.  Parabolic  Quarter  Thickness  of  Top 
of  K2  Region  vs  Height  of  Peak  of  F2  Region, 
Individual  observations  from  incoherent  scatter  at 
Millstone  Hill;  nighttime 
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Figure  2c.  Parabolic  Quarter  Thickness  of  Top 
of  F2  Region  vs  Height  of  Peak  of  F2  Region. 
Individual  observations  from  incoherent  scatter  at 
Millstone  Hill:  near  sunrise  anci  sunset 
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145) 
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max 


is  a simple  engineering  approximation.  The  choice  of  145  km  as  reference  point 
is  arbitrary  but  consistent  wit'n  both  Figures  1 and  2. 

To  develop  an  empirical  model  around  such  a relationship,  all  the  data  were 
pooled  and  then  the  mean  k of  Eq.  (1)  was  determined  for  cells  of  size  month  by 
hour.  Results  are  contoured  in  Figure  3. 


Z\WH  2IJUN  23  SEP  21  DEC 


By  inspection,  the  diurnal  variation  and  the  seasonal  variation  of  k are  not  sep- 
arable; consequently,  the  seasonal  mean  and  the  diurnal  mean  of  the  total  sample 

7 

are  of  no  direct  assistance  in  modeling  the  parameter.  Becker  has  already  shown 
a dependence  c f on  the  mean  solar  flux  (^)  and  the  noon  solar  zenith  angle  (x>: 


or  cos  X 


Unfortunately,  the  number  of  data  points  per  month  per  year  of  this  set  is  small, 
so  in  this  paper  ♦ is  assumed  constant. 

7.  Becker,  W.  (1969)  The  seasonal  anomaly  of  the  F region  at  mid  latitudes  and 

its  interpretation,  in  Electron  Density  Profiles  in  Ionosphere  and  Exosphere, 
Jon  Frihagen,  Ed.  , North  Holland  Publishing  Company,  Amsterdam,  'ils-'iSO. 
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With  this  insii^ht.  a trial  model  discontinuous  at  sunrise  (SH)  and  sunset  (SS) 
was  constructed  such  that  the  value  of  k in  tlq.  (1)  was 

k - 0.  34 

for  the  night  hours  (SS  < t < SR),  or  (2) 

k = 0.  25  ( 1 + cos  x) 
for  the  day  hours  (SR  ^ t s SS), 

where  t is  local  time.  The  difference  between  modeled  k and  the  observed  mean  k 
is  contoured  in  F.lgure  4.  Background  errors  are  generally  less  than  10  percent, 
while  the  worst  errors  are  about  25  percent  in  summer  and  about  20  pei  cent  in 
winter. 

aiMAR  2IJUN  23  SEP  21  DEC 


Figure  4,  Difference  Between  Modeled  k and  Observed  k 
for  Eq.  (2) 

If  empirical  terms  are  added  for  the  systematic  variations  in  the  transition 
f>eriods  near  sunrise  and  sunset,  then  a relation  such  as 


10 


k 


0.045 


<3) 


^ 0.090  cos  X , ) 0. 34  ifSS<t<SR 

(1  + (t-SS-1)^)  (1  + (t -SR -0.7)^  I 0.25(1  + cos  x>  ifSR^t^SS 

273 

can  model  mean  k with  residual  errors  as  shown  in  Figure  5. 


Figure  5.  Residual  Errors  in  Modeling  k When  Simple  Terms  are 
Added  for  Sunrise  and  Sunset  as  in  Eq  (3) 


Inasmuch  as  it  is  proposed  that  this  simple  model  be  applied  to  other  stations 
at  higher  and  lower  latitudes,  as  suggested  by  Figiire  1,  it  would  not  seem  proper 
to  model  the  prominent  features  found  during  nighttime  in  winter  and  during  the 
early  day  in  summer,  since  they  may  well  be  local  effects. 


3.  TEST 

An  estimate  of  the  usefulness  of  the  proposed  model  is  made  by  testing  it 
against  the  individual  observations  from  which  it  was  derived.  As  a test  parameter, 
we  use  variability  defined  as  the  sum  of  the  squares  of  the  deviations  between  the 
observed  values  and  predicted  values. 


’1 


For  a first  test,  assume  that  some  external  climatology  could  provide,  with- 
out bias,  estimates  of  the  mean  value  of  the  topside  thickness  for  each  of  the 
samples  of  Fig\u*e  2.  Then  Table  1 shows  that  the  specification  model  suggested 
here,  using  an  accurate  value  of  can  provide  individual  estimates  of  topside 

thickness  during  the  daytime  with  about  56  percent  less  variability  than  the  varia- 
bility about  the  sample  mean.  Note,  however,  that  while  the  central  value  of  the 
night  sample  has  been  specified  usefully,  the  technique  has  not  reduced  the  spread 
of  individual  values. 

Table  1.  Variability  for  Samples  of  Figure  2 


Simrise  to 
Sunset 

Sunset  to 
Sunrise 

Sunrise  t 2 h 
Simset  ± 2 h 

Sample  size 

1795 

1608 

1141 

Avei^ge  variability  about 
sample  mean 

256 

208 

263 

Average  variability  using 
n.odel 

112 

187 

134 

Reduction  of  variability 
using  model 

56% 

10% 

49% 

Consider  a climatology  that  could  provide  mean  values  as  a function  of  both 
month  and  time  of  day;  then  the  variability  might  not  be  reduced  as  much  as  sug- 
gested by  Table  1.  To  test  against  such  a climatology,  we  divided  the  observations 
into  cells  2 months  by  2 hours  (about  50  data  points  per  cell).  Linear  regression 
relations  for  each  cell  were  compirted  to  approximate  an  efficient  specification 
method.  Results  are  shown  in  Table  2,  The  variability  has  been  computed  with 
reference  to  the  mean  value  of  each  cell.  During  winter  (Nov.  Dec.  Jan.  Feb)  and 
equinox  (Mar.  Apr.  Sep.  Oct)  days,  the  simple  model  presented  here  is  useful  and 
not  significantly  inferior  to  the  efficient  relations  derived  from  linear  regression. 
As  expected,  the  model  fails  significantly  for  summer  mornings  and  winter  nights, 
those  periods  where  k of  Eq.  (3)  was  known  to  differ  from  the  observations.  It  is 
of  general  interest  to  add  that  the  low  values  of  reduction  of  variability  shown  by 
the  linear  regression  results  for  these  same  periods  implies  that  the  physical 
mech^htism  which  correlates  individual  observations  of  stfid  is  of  negligible 

importance  during  these  periods. 

The  technique  should  be  verified  with  an  independent  source  of  data. 


12 


0) 


a 

E 


ctf 

03 


O 

.c 


X 

'S 

o 

E 

ca 


u 

o 


u 

> 


c 

o 


<J 

SJ 

-o 

43 


4) 

•§ 

H 


beg 

3 

< 

AA  A 

mcsoirjMCTivcoMcoc'ieo 

O'^'-iooo'^CMiooonc'i 

• • • vn/  * * * v * 

>-3  P5 

0) 

iA0)t*(0C4O00C0tAl^00C^l 

P 

•-3 

1 

AAAAA  A A 

0)OC^OO(D^OlO«OCOO> 

Scs 

^e0^<00<0(000)t^^€<) 

CMOOOOO^COOOO^ 

o 

o 

AA  AA  AAA 

000)10>-<m<OMCMCO(0^0 

OOCMOO’-'M'^'^VCJ^ 

V/V  *W  * * • *VW 

c 

0 

T 

a 

43 

§ 

03 

W 

■3 

O' 

M 

1= 

A A 

<Oi>.<ooQOa)Oooc4t<*i>c^) 

V V * 

<3 

Sofi 

<^COO^<COOO>lOO>COCO 

WO^^’-^OJ^CMCOCOO’-* 

43 

1 

C 

»-3  CP 
u 

AAA  A A 

C'jirt^r-eotrtC5«o^^icco 

^»-<ioeoc^c4wco^^»-^0 

vvn; V-* 

0»-»000»-‘<000t^^^?0'^ 

■£ 

V S 

43 

Q 

A A A AA 

ooo^^r-eco^»ot^05 

OOCO'^^^COCOO’-'OO 

V V V vV 

L 

^<ocofO<ocaoo»-'<o^ 

ooo^ift^co^cacg^>o 

Local 

Standard 

Time 

C4^<OOOOC'l^<OCCOW'^ 

1 1 1 1 1 I 1 1 1 • 1 1 

OC^^COOOOC4^<OOOOC^ 

13 


* First  entry  percentage  reduction  of  the  variability  using  linear  regression  (R)  compared  to  sample  mean; 
second  entry  percentage  reduction  or  increase  of  variability  using  model  (M)  compared  to  sample  mean. 
An  increase  is  denoted by  carets  «>). 


4.  DISCUSSION 


From  either  climatology  or  a local  measurement  h is  easily  obtained;  a 
thickness  parameter  is  not.  The  useful  correlation  displayed  between  the  two 
parameters  over  day,  season,  solar  cycle  and  latitude  recommends  the  use  of 
a simple  specification  technique  such  as  the  one  presented  here.  There  are  indi- 
cations in  both  Figures  1 and  2 that  a second  order  relation  in  h „ would  be  more 
exact,  but  unnecessary  at  this  stage.  In  fact,  to  first  order.  Eq.  (2)  is  sufficient 
since  the  improvement  using  the  add’tional  terms  of  Eq.  (3)  is  only  about  10  per- 
cent. 

While  this  technique  may  be  useful  in  modeling  the  upper  ionosphere  for  vari- 
ous engineering  applications,  the  significance  of  the  basic  physics  contained  in 
Eq.  (1)  is  far  from  clear.  The  absence  of  appreciable  diurnal  or  seasonal  varia- 
tion of  the  thickness  of  the  lower  E region  is  generally  accepted,  but  note  that  the 
choice  of  145  km  as  the  null  height  of  Eq.  (1)  was  arbitrary  and  would  be  changed 
if  a second  order  relation  were  substituted.  On  the  other  hand,  it  seems  to  be 
consistent  with  many  observations  that  suggest  the  stability  of  atmospheric  param- 
eters in  the  vicinity  of  about  120  km. 

Previous  analysis  of  various  models  of  the  upper  atmosphere  suggests  that  the 
thickness  parameter  varies  directly  wi-.h  the  temperature  of  the  electrons  and  ions 
and  inversely  with  the  force  of  gravity  and  the  mean  molecular  mass.  The  varia- 
tion of  the  force  of  gravity  with  height  is  a very  small  portion  of  the  variation 
suggested  by  Figure  2;  direct  observations  by  satellite  probes  suggest  that  at  the 
heights  of  interest,  between  200  and  500  km,  the  mean  molecular  mass  is  very 

4- 

nearly  that  of  0 . Therefore,  the  variation  of  the  thickness  of  either  the  upper  or 
lower  F region  is  an  expression  of  the  effective  ion-electron  temperature.  From 
this  it  is  concluded  that  either  the  thickness  or  the  height  of  the  F region  is  a 
macroparainetric  measure  of  the  temperature. 

The  change  of  thickness  parameter  with  the  day  to  day  change  of  solar  or  geo- 
physical activity  is  considered  a promising  separate  study. 
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